function plot_calibrated_bellows_toolbox(patient, fig_num)
    
     %| Plot original vs drifted bellows voltage, using the patient.drift
     %| as;
     %|     patient.bellows_volt_drifted = patient.bellows_volt + patient.drift*patient.scan_times_zerod; 
    
    %------------------------------------------------------------------------
    %|      Dependancies;  
    %|                  
    %|                 
    %|                 
    %------------------------------------------------------------------------
    %   This file is part of the
    %   5D-Novel4DCT Toolbox  ("Novel4DCT-Toolbox")
    %   DH Thomas, Ph.D
    %   University of California, Los Angeles
    %   Contact: mailto:dhthomas@mednet.ucla.edu
    %------------------------------------------------------------------------
    % $Author: DHThomas $	$Date: 2014/04/01 10:23:59 $	$Revision: 0.2 $
    
    if nargin < 2
        figurecount = length(findobj('Type','figure'));
        fig_num = figurecount+1;
      
    end
    
    bellows_calibrated = polyval(patient.volume_calibration, patient.bellows_volt_drifted);
    flow_calibrated = polyval(patient.volume_calibration, patient.flow_drifted);
    
    volume_5pctile=prctile(bellows_calibrated(:),5);
    volume_25pctile=prctile(bellows_calibrated(:),25);
    volume_50pctile=prctile(bellows_calibrated(:),50);
    volume_75pctile=prctile(bellows_calibrated(:),75);
     volume_95pctile=prctile(bellows_calibrated(:),95);
   

    figure(fig_num);
        set(gcf,'units','normalized','position',patient.default_figure_position);

    subplot(2,1,1)
    hold on;
    p1 = plot(patient.scan_times_zerod(:),bellows_calibrated(:),'.');
    p5 = plot([0 patient.scan_times_zerod(end)],[volume_5pctile volume_5pctile],'r--', 'LineWidth',3);
    p25 = plot([0 patient.scan_times_zerod(end)],[volume_25pctile volume_25pctile],'g--','LineWidth',3);
        p50 = plot([0 patient.scan_times_zerod(end)],[volume_50pctile volume_50pctile],'r--','LineWidth',3);
    p75 = plot([0 patient.scan_times_zerod(end)],[volume_75pctile volume_75pctile],'g--','LineWidth',3);
    p95 = plot([0 patient.scan_times_zerod(end)],[volume_95pctile volume_95pctile],'r--','LineWidth',3);

    legend([p5(1,1);p95(1,1)], '5th %ile', '95th %ile')
    hold off
    title('Breaths sampled during scans')
    xlim([0 patient.scan_times_zerod(end)])
    xlabel('Scan Time [s]')
    ylabel('Lung Volume [L]')
    
        subplot(2,1,2)
    hold on;
    p1 = plot(bellows_calibrated,'b');
    p5 = plot([0 patient.dim(3)],[volume_5pctile volume_5pctile],'r--', 'LineWidth',3);
        p5 = plot([0 patient.dim(3)],[volume_5pctile volume_5pctile],'r--', 'LineWidth',3);
    p25 = plot([0 patient.dim(3)],[volume_25pctile volume_25pctile],'g--','LineWidth',3);
        p50 = plot([0 patient.dim(3)],[volume_50pctile volume_50pctile],'r--','LineWidth',3);
    p75 = plot([0 patient.dim(3)],[volume_75pctile volume_75pctile],'g--','LineWidth',3);
    p95 = plot([0 patient.dim(3)],[volume_95pctile volume_95pctile],'r--','LineWidth',3);

    legend([p5(1,1);p95(1,1)], '5th %ile', '95th %ile')
    hold off
    title('Breaths sampled during scans')
    xlim([0 patient.dim(3)])
    xlabel('Slice Number [mm]')
    ylabel('Lung Volume [L]')
